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Abstract For star-forming regions, there is a correlation of radio and FIR-emission 
established. The radio emission is caused by synchrotron radiation of electrons, 
while the FIR emission is attributed to HII regions of OB stars and hot dust powered 
by those stars. Another possible correlation for star-forming regions might exist be- 
tween the aforementioned radio emission and the gamma ray emission in the GeV 
regime. The GeV gamma ray emission of star-forming regions is dominated by the 
decay of neutral pions formed in collisions of cosmic ray (CR) protons accelerated 
at supernova remnants (SNRs) with ambient hydrogen, while the electrons gener- 
ating the synchrotron emission are assumed to be accelerated at the same SNRs. 
Assuming the same spectral shape for the CR proton and electron distribution and 
a fixed ratio of CR protons to electrons, the proton- and electron-calorimetry of the 
star-forming regions can be tested. Furthermore, this provides a method to derive 
the magnetic field strength in the star-forming region complementary to the existing 
methods. 



1 Introduction 



A plot of the total gamma ray luminosity against the luminosity at 5 GHz for 
M82 and NGC 253 (both starburst galaxies), NGC 1068 and NGC 4945 (starburst- 
Seyfert composites) and 30 Doradus in the Large Magellanic Cloud and the Milky 
Way center reveals a trend, as Fig. [T] shows. To interpret this possible correlation, 
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Fig. 1 Correlation between radio luminosity at 5 GHz and the y-luminosity. 



a theoretical approach is followed: The emission of the objects dealt with is as- 
sumed to be dominated by SNRs inside the objects. The detected radio emission is 
attributed to synchrotron radiation of electrons accelerated at the SNRs, the GeV 
gamma rays are caused by the decay of neutral pions formed in collisions of protons 
accelerated at the SNRs with ambient protons: 



P CR+p^7T +/ - /0 +X (1) 

7i° -^y+y. (2) 

Neutrinos are produced from charged pions, where searches are dominated by a 
background of atmospheric neutrinos so far [2|. Both protons and electrons are 
probably accelerated at the same site, so it is assumed that their spectra have the 
same -or a very similar- shape, they only differ by a constant factor 



^=J E A' P ,^=4^)", 0.01 <**<!. (3) 
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To examine the trend in Fig. [T| both radiation processes are investigated in more 
detail in the following sections. 
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2 Electron synchrotron radiation — > Radio emission 
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The electron spectral energy distribution is assumed to be a simple power-law: 

N{E e ) dE e = a e (-^j) dE . (4) 
\m e c z J 

If each electron radiates all its energy at a single frequency v, the emission coeffi- 
cient £ v can be written as: 



E = Y m e c 2 , 7= f — J , — = ^-m e c 2 , % = — , 8 = {\-p)/2, (6) 

2v G 

using the expressions above, where e is the elementary charge, Vg is the gyrofre- 
quency and at is the Thomson cross-section. In order to obtain the radio luminosity 
of an object, the emission coefficient is multiplied by the average volume V of an 
SNR and the number of SNRs, Nsnr in the object: 

V] 

L Radio = vn snr J e v (v)dv (7) 

vo 

dL v =VN M \f> 2 my (J-Y dv . (8) 

By expressing the emission coefficient in terms of energy instead of frequency, the 
synchrotron flux can be written as: 

*£-!^-t^(?)(£)'- "> 

where A e = VWsNRfle- 



3 GeV gamma rays from pp-interaction 

Following [4], the differential proton flux in each object is derived from gamma ray 
observations under the assumption that these gamma rays are dominantly formed 
by pion decay from proton-proton interactions. The spectral shape of the primary 
proton flux is the same as that of the primary electrons (see eq.[T]i and of the gamma 
rays from pion decay in the energy range of 300 MeV < Ey < E^ ax /10. The 



4 



S. Schoneberg, J. Becker Tjus, F. Schuppan 



differential gamma ray luminosity is calculated as: 

= VN SNR a p n H cE^ J a lnel (£ p ) (^j ' F y ^ • (10) 

^jnin 

The calculation is done following |4|, where the 5 -functional approach is used for 
energies below 100 GeV and the analytical approximation given there is used for 
higher energies. 



4 Comparison and conclusions 

Combining the spectra as derived from protons and electrons, respectively, allows 
to model the SED for each object examined. In Fig. |4] this is shown for M82. The 
gamma ray luminosity Ly and the spectral index F for NGC 4945, 1068 and 3043 
given by and for NGC 253 from JT| are used to calculate A p = VA^sNRflp- Using 
radio data, fits to the radio spectra are made to find the radio synchrotron index 8 
and the normalization of the primary electron spectrum, A e , for a magnetic field 
of B = 1 jiG. From 8, the spectral index p of the underlying primary electron 
spectrum is derived. The corresponding values of the calculation are given in Table 
H where «h = 100 cm~ 3 was assumed for the calculation of A p . The quantity A e 
depends on the magnetic field as A e °< B s /B, whereas A p depends on the average 
hydrogen density «h as A e °= 1 /«h- These expressions can be combined using the 
inequality in eq. (3). This allows to determine a possible range of combinations of 
the magnetic field strength B and the average hydrogen density «h- This is done in 
Fig. [2] As can be seen, for each galaxy the observed values for the average hydrogen 
density and the magnetic field strength lie within the theoretically allowed range, 
which is the area between the lines corresponding to a e /a v = 0.01 (dashed black 
line) and a e /a p = 1 (solid red line). To test the assumption that the energy spectra of 
protons and electrons have the same shape, the electron and proton spectral indices 
of the individual galaxies have been calculated and are shown in Fig. [3] For a given 



Source 


r=p v 


A p [l/ergj 


8 


Pe 


A e [l/erg] 


NGC 1068 
NGC 4945 
NGC 253 
NGC 3034 


2.31±0.13 
2.3±0.1 

2.14±0.18 
2.2±0.2 


1.45E+68 
6.16E+67 
l.OOE+66 
3.20E+66 


-0.710± 0.004 
-0.564± 0.011 
-0.525 ±0.006 
-0.457 ±0.001 


2.421 ±0.008 
2. 129 ±0.022 
2.051 ±0.012 
1.914 ±0.002 


1.76E+69 
7.50E+66 
1.77E+66 
9.60E+65 



Table 1 Parameters A p calculated from F-index and gamma ray luminosity, A e calculated from a 
fit to the radio data. 



average hydrogen density and an estimate of a e /a p , this approach could be used 
to determine the magnetic field strength, complementary to the methods presented 
in IS] and 0. For a given average hydrogen density, magnetic field strength and 
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Fig. 2 Possible range of the parameters B and «h for the individual galaxies. Density estimates 
taken from: f7l[6). 



a fixed ratio for a e /a p , this approach could be used to check whether the galaxy 
is calorimetric for both protons and electrons to the same extent. However, further 
improvement of the statistics is required to check the possible correlation of GeV 
gamma rays and radio emission from star-forming regions. Additionally, a more 
detailed theoretical examination has to be performed to establish this connection 
and gain new insights. The versatility of such a correlation strongly motivates further 
investigation in the future. 
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Fig. 3 Spectral index of the proton and electron distribution for each object as derived from gamma 
or radio data, respectively. The solid line represents p p = p e ■ 
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Fig. 4 Modeled partial SED of M82. 
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